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Abstract 

We consider the Cauchy problem for the massless scalar wave equation in the Kerr 
geometry for smooth initial data compactly supported outside the event horizon. We 
prove that the solutions decay in time in Lfjf c . The proof is based on a representation 
of the solution as an infinite sum over the angular momentum modes, each of which 
^vO ' is an integral of the energy variable ui on the real line. This integral representation 

involves solutions of the radial and angular ODEs which arise in the separation of 
1^ variables. 
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1 Introduction 



In this paper we study the long-time dynamics of massless scalar waves in the Kerr ge- 
ometry. We prove that solutions of the Cauchy problem with smooth initial data which is 
compactly supported outside the event horizon, decay in L™ c . Our starting point is the 
integral representation for the propagator [5], which involves an integral over a complex 
contour in the energy variable u>. In order to study the long-time dynamics, we must 
deform the contour to the real line. To this end, we carefully analyze the solutions of 
the associated radial and angular ODEs which arise in the separation of variables. In 
particular, we show that the integrand in our representation has no poles on the real axis. 
We call such poles radiant modes, because in a dynamical situation they would lead to 
continuous radiation coming out of the ergosphere. 

We now set up some notation and state our main result. As in [5], we choose Boyer- 
Lindquist coordinates (t, r, tp) with r>0,0<i9<7r,0<<y9<27r, in which the Kerr 
metric takes the form 

ds 2 = A (rft _ as . n 2# d(p) 2 _ U f^ + d ^\ _ *^( a dt _ (r 2 + a 2 )d^) 2 (1.1) 
with 

U{r,<&) = r 2 + a 2 cos 2 tf, A(r) = r 2 - 2Mr + a 2 , 

where M and aM denote the mass and the angular momentum of the black hole, respec- 
tively. We restrict attention to the non-extreme case M 2 > a 2 , where the function A has 
two distinct zeros, 



r = M - \jM 2 - a 2 and n = M + V ' M 2 - a 2 , 

corresponding to the Cauchy and the event horizon, respectively. We consider only the 
region r > n outside the event horizon, and thus A > 0. The ergosphere is the region 
where the Killing vector is space-like, that is where 

r 2 -2Mr + a 2 cos 2 $ < 0. (1.2) 

The ergosphere lies outside the event horizon r = r\ , and its boundary intersects the event 
horizon at the poles $ = 0, 7r. 

Theorem 1.1 Consider the Cauchy problem for the wave equation in the Kerr geometry 
for smooth initial data which is compactly supported outside the event horizon and has 
fixed angular momentum in the direction of the rotation axis of the black hole, i.e. for 
some k £ Z, 

($o,3t$o) = e-** (*o,$*o)M) G Co°°((Vi, oo) x S 2 ) 2 . 



Then the solution decays in L?° J(ri, oo) x S ) as t 



oo. 



The study of linear hyperbolic equations in a black hole geometry has a long history. 
Regge and Wheeler [12] considered the radial equation for metric perturbations of the 
Schwarzschild metric. In the late 1960s and early 1970s, Carter, Teukolsky and Chan- 
drasekhar discovered that the equations describing scalar, Dirac, Maxwell and linearized 
gravitational fields in the Kerr geometry are separable into ordinary differential equations 
(see [2]). Much research has been done concerning the long-time behavior of the solutions 
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of these equations, through both numerical and analytical methods. Price [TU] gave argu- 
ments which indicated decay of solutions of the scalar wave equation in the Schwarzschild 
geometry. Press and Teukolsky [9j did a numerical study which strongly suggested the 
absence of unstable modes, and Whiting [IT] later proved that for uj in the complex plane, 
such unstable modes cannot exist. This "mode stability" does not rule out that there might 
be unstable modes for real uj (what we call radiant modes). Furthermore, mode stability 
does not lead to any statement on the Cauchy problem. Finally, Kay and Wald [7j used 
energy estimates to prove a boundedness result for solutions of the scalar wave equation 
in the Schwarzschild geometry. 

Unfortunately, these energy methods cannot be used in a rotating black hole geometry, 
because the energy density is indefinite inside the ergosphere, making it impossible to 
introduce a positive definite conserved scalar product. This difficulty was dealt with in [H 
[6], where Whiting's mode stability result was combined with estimates for the resolvent 
and for the radial and angular ODEs. In [5] we established an integral representation 
which expresses the solution as a contour integral of an integrand involving the separated 
radial and angular eigenfunctions over a contour staying within a neighborhood located 
arbitrarily close to the real axis. This integral representation is the starting point of the 
present paper. After deforming the contours onto the real axis, we can prove decay using 
the Riemann-Lebesgue Lemma, similar to the case of the Dirac equation [1]. We remark 
that the decay result of our paper would not be expected to hold for a massive scalar field 
satisfying the Klein-Gordon equation, as indicated in [TJ. 

Finally, we note that the problem considered here is closely related to one of the 
major open questions in general relativity; namely the problem of linearized stability of 
the Kerr metric. For the stability under metric perturbations one considers the equation 
for linearized gravitational waves, which can be identified with the general wave equation 
for spin s = 2 (see [2]). Thus replacing scalar waves (s = 0) by gravitational waves 
(s = 2), the above theorem would prove linearized stability of the Kerr metric. However, 
the analysis for s = 2 would be considerably more difficult due to the complexity of the 
linearized Einstein equations. Nevertheless, we regard this paper as a first step towards 
proving linearized stability of the Kerr metric. 



2 Preliminaries 

We recall a few constructions and results from [5] [6] which will be needed later on. As 
radial variable we usually work with the Regge- Wheeler variable defined by 

du _ r 2 +a 2 

dr ~ A ' 1 j 

then u = — oo corresponds to the event horizon. It is most convenient to write the wave 
equation in the Hamiltonian form 

id t y = H^, (2.2) 
where ^/ = (<&,idt&). The Hamiltonian can be written as 

h= (°aI)- (23) 
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where 



A 

(3 



8 2 ^ 2 5 

<9u 9u 



A 



2i.2 n 



2a/c 
P 

2 , 2 

r + a 



1 



A 



r 2 + a 2 



a 2 sin 2 1? 



-jT - ! — 2 A S 2 

r z _j_ a z 



A 



a 



r 2 + a 2 



r 2 + a 2 



(2.4) 
(2.5) 
(2.6) 



The operators A and j5 are symmetric on the Hilbert space L 2 (M. x S 2 ,dfi) 2 with the 
measure 

<i/i := pdud cos??. (2-7) 



It is immediately verified that the Hamiltonian is symmetric with respect to the bilinear 
form 

T T f /T ( A 

JRxS 2 



1 



* 2 }c 2 dji . 



[2.1 



As is worked out in detail in [5], the inner product <\I',\E'> is the physical energy of ^ . 
Therefore, we refer to <., .> as the energy scalar product. The fact that the energy scalar 
product is not positive definite can be understood from the fact that the operator A is 
not positive on L 2 (R x S 2 ,d/i). 
Using the ansatz 

(2.9) 



$(t, r, i?, (p) = e -*>*-*l«P R( r ) 0(0), 



the wave equation can be separated into an angular and a radial ODE, 



(2.10) 



lo £ C I lima; I < 



Here the angular operator A<j fe is also called the spheroidal wave operator. The separation 
constant A is an eigenvalue of »4 W) fc and can thus be regarded as an angular quantum 
number. In |6j it was shown that if a; is in a small neighborhood of the real line, more 
precisely if 

1 + |Rew| J ' 

then for sufficiently small e > the angular operator A^^ has a purely discrete spec- 
trum (A n ) ng N with corresponding one-dimensional eigenspaces which span the Hilbert 
space L 2 (S 2 ). We denote the projections onto the eigenspaces by Q n (k,uj). These projec- 
tions as well as the corresponding eigenvalues A n are holomorphic in u £ U £ . In analogy 
to the eigenvalues 1(1 + 1) of the Laplacian on the sphere, the angular eigenvalues A n grow 
quadratically for large n in the sense that there is a constant C(k,u) > such that 



\X n (k,u)\ > 



n 



C(k,u) 



for all n G N 



We set 



ak 



uJ() 



ri + a 1 



with r\ the event horizon and use the notation 

£1(uj) = lu — LOQ . 



(2.11) 



(2.12) 



(2.13) 
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In order to bring the radial equation into a convenient form, we introduce a new radial 
function <p(r) by 

(j)(r) = V ' r 2 + a 2 R(r) . 

Then in the Regge- Wheeler variable, the radial equation can be written as the "Schrodinger 
equation" 

-* + V(u))(t>(u) = (2.14) 



du 2 

with the potential 

V M = -L + -^X + ^^ + T^^V 7 ^. (2-15) 
V r 2 + a 2 J (r 2 + a 2 ) 2 ^ r 2 + a 2 

In [5] we derived an integral representation for the solution of the Cauchy problem of 
the following form, 



27rz^ e ^ e\0 \Jr Jc- 

k£Z neIN VJUe J ° e 



duj e~ iwt (Q k , n (u) Soo(lo) tfg)(r,0) . (2.16) 
Here the integration contour C £ must lie inside the set U £ . 



3 Asymptotic Estimates for the Radial Equation 
3.1 Holomorphic Families of Radial Solutions 

In this section we fix the angular quantum numbers k, n and consider solutions (f> and </> 
of the Schrodinger equation (I2.14|) which satisfy the following asymptotic boundary con- 
ditions on the event horizon and at infinity, respectively, 

lim e~ inv 4>(u) = 1, lim ( e~ inu 4>(u))' = (3.1) 

lime iuJU 0(u) = 1, lim (e^^u))' = 0. (3.2) 

These solutions were introduced in [5j for uj in the lower complex half plane intersected 
with U e . Here we will show that they are holomorphic in uj, and we will extend their 
definition to a larger w-domain. More precisely, we prove the following two theorems. 

Theorem 3.1 The solutions 4> are well-defined on the domain 

D = U £ n i lo G C I Imu < ' ' ~ "' 



2(rf + a 2 

They form a holomorphic family of solutions in the sense that for every fixed u G M 
and n£N, the function 4>(u) is holomorphic in uj G D. 

Theorem 3.2 For every angular momentum number n there is an open set E containing 
the real line except for the origin, 

E D E := U £ n {uj G C | Imu < and uj ^ 0} , (3.3) 

such that the solutions <j) are well-defined for all uj G E and form a holomorphic family 
on E. 
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For the proofs we will rewrite the Schrodinger equation with boundary conditions (|3.1l 
13. 2p as an integral equation (which in different contexts is called Lipman-Schwinger or 
Jost equation). Then we will perform a perturbation expansion and get estimates for all 
the terms of the expansion. To introduce the method, we begin with the solutions <j>; the 
solutions (j) will be treated later with a similar technique. First we write the Schrodinger 
equation (|2.14p in the form 

- ku) = -W{u) (3.4) 

with a potential W = + V(u) which vanishes at u = —oo. We define the Green's 
function of the differential operator — d 2 — VL 2 by the distributional equation 

(-d 2 - CI 2 ) S(u, v) = 5(u-v). (3.5) 

The Green's function is not unique; we choose it such that its support is contained in the 
region v < u; i.e. 

S(u,v) = @(u-v) x < 2ifl\ J (3.6) 

[ v-u if n = . 

(Here denotes the Heaviside function defined by G(x) = 1 if x > and @(x) = 
otherwise.) We multiply (|3.4ft by the Green's function and integrate, 

S(u, v) ((-d 2 - n 2 ) (0(u) - e inu )) dv = - ! S(u, v) W{v) (f>(v) dv . 

J — oo 

If we assume for the moment that (j) satisfies the desired boundary conditions (|3.1|) . we can 
integrate by parts on the left and use (|3.5[) . This gives the Lipman-Schwinger equation 

/u 
S(u,v)W(v)cj>(v) dv, 
-oo 

which in the context of potential scattering is also called Jost equation (see e.g. [3]). Its 
significance lies in the fact that we can now easily perform a perturbation expansion in 
the potential W. Namely, taking for <j> the ansatz as the perturbation series 



(3.7) 



1=0 

we are led to the iteration scheme 

0(°) (u) = e inu 
^C+ 1 )(u) = - / S(u,v) W(v) k l \v) dv 



(3i 



This iteration scheme can be used for constructing solutions of the Jost equation, and this 
will give us the functions <p with the desired properties. 

Proof of Theorem I3.il Fix to £ D. As the potential W is smooth in r and vanishes on 
the event horizon, we know that W has near n the asymptotics W = 0(r — n). This 
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means in the Regge- Wheeler variable (|2.ip that W decays exponentially as u — ► — oo. 
More precisely, there is a constant c > such that 

\W(u)\ < ce< u with 7 := r \ - r ° . (3.9) 

rf + a 2 

Let us show inductively that 

< ,.'«--- with „:. (7 _ Im „_ |Imii|)2 . (3.10) 

In the case Z = 0, the claim is obvious from (13. 8h , Thus assume that A3. lQf) holds for a 
given /. Then, estimating the integral equation in (|3.8p using (|3.9p . we obtain 

< cfjf [ U \S(u,v)\e^- lmn)v dv . (3.11) 

J — oo 

The Green's function (|3.6p can be estimated in the case u < u by 



e -iO(«-lO r ^ 



< ( u - u ) e l Imn l( u -») 



2 

Substituting this inequality in (13.111) gives 

I j, (7+1) / M ^ / llmniu \ ( \ (7— Im f2— llm Sll) u j 

'[u)\ < c 11 e 1 / [u — v) e w 1 u av . 



Since the parameter a := 7 — Imfi — |Im$7| is positive according to the definition of D, 
we can carry out the integral as follows, 



d 

da J ./-oo V da J a a 2 



(u - v) e av dv = ( u - — ) j e av dv = [u ; , 



This gives ()3.10j) with / replaced by I + 1. 

Since for u on a compact interval, the analytic dependence of the solutions in u from 
the coefficients and the initial conditions follows immediately from the Picard-Lindelof 
Theorem, it suffices to consider the region u < uq for any uq £ M. By choosing no suffi- 
ciently small, we can arrange that u < 1/2 for all u < uq. Then the estimate (|3.10j) shows 
that the perturbation series (13. 7p converges absolutely, uniformly inuS (— 00, uo). Using 
similar estimates for the n-derivatives of <j)( 1 ' , one sees furthermore that the perturbation 
series (I3.10P can be differentiated term by term, and using (I3.5P we find that (ft is indeed 
a solution of (|3.4p . Furthermore, 



<P(u)-e^ u = X> (/) (u), 



1=1 

and taking the limit u — > 00 and using (|3.10p we find that the right side goes to zero. 
Using the same argument for the first derivatives, we obtain (|3.ip . 

In order to prove that (ft is analytic in u, we first note that if 0, 7^ 0, we can differ- 
entiate the perturbation series (13. 7p term by term and verify that the Cauchy-Riemann 
equations are satisfied (note that A n is holomorphic in to according to |6j). Since <ft is 
bounded near = 0, it is also analytic at O = 0. ■ 
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We turn to the solutions (ft. In analogy to (|3,4p . we now write the Schrodinger equation 

as 2 

-uA <Mu) = -W(u)j>{u) (3.12) 
with 

, . 2ak (ah) 2 An A 1 ~ 

W(u) = -u^—^ - -, l' 2 , 2 + ? " 2 . 2 + <9 2 vW a 2 . 3.13 

Assuming that u; ^ 0, we choose the Green's function as 

S(u,v) = — ( e' iuJ ^- u) - e iuj{v ' u) ) Q(v - u) . (3.14) 

The corresponding Jost equation is 

;>oo 

^(n) = e -*™- / S(u,v)W(v)4>(v)dv. 

J u 

The perturbation series ansatz 



E^ (0 ( 3 - 15 ) 



z=o 



leads to the iteration scheme 



(3.16) 

Note that, in contrast to the exponential decay (|3.9p . now the potential W-^, (|3.13|) . has 
only polynomial decay. As a consequence, the iteration scheme allows us to construct <fi 
only inside the set Eq as defined in 



Lemma 3.3 The solutions (j) are well-defined for every u E Eq. They form a holomorphic 
family in the interior of Eq. 

Proof. Fix to € Eq. Then u ^ and Imw < 0, and this allows us to estimate the 
potential (|3.13|) and the Green's function (|3.14p for u, v > uq and some uq > by 



\W(v)\ < 4, \S(u,v)\ < J-e 1 ™"^ . (3.17) 

v z \u>\ 

Let us show by induction that 

\^ l \u)\ < If*)*,**" «. 
l\ \\uj\uj 

For I = this is obvious from (|3.16|) . whereas the induction step follows by estimating the 
integral equation in (|3.16p with (|3.17p . 



^ (,+1)wl 5 Km) f h 



Imcj (u— v) Imijjv 
v 2+l e 



(I + 1)! V>| it 
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Hence the perturbation series (|3.15p converges absolutely, locally uniformly in u. It is 
straightforward to check that <p satisfies the Schrodinger equation (|3,12p with the correct 
boundary values f|3.2[) . If Imu < 0, one can differentiate the series (|3. 15|) term by term 
with respect to uj and verify that the Cauchy-Riemann equations are satisfied. ■ 



It remains to analytically extend the solutions <j) for fixed n to a neighborhood of any 
point ujq E R \ {0}. To this end, we need good estimates of the derivatives of <j> with 
respect to uj and u. It is most convenient to work with the functions 

^ l \u) := (2ioo) 1 e iuJu ^ l \u) , (3.18) 

for which the iteration scheme (|3.16p can be written as V> (0) = 1 and 

^ +1 )(u) = / (e-^M -1) W(v)i/>®(v)dv. (3.19) 

J u 



Lemma 3.4 For every wo £ 1\ {0} and n E N ; there are positive constants c, K, 5, such 
that for all uj E Eq n B$(ujq) with Imw < and all p,q,n E N the following inequality 
holds, 

d V ( d 



(3.20) 



Proof. According to [6j, A n is holomorphic in a neighborhood of ujq, and thus (for example 
using the Cauchy integral formula) its derivatives can be bounded in B$(ujq) by 



for suitable K > 0. Since the potential W, (|3.13p . is also holomorphic in r (in a suitable 
neighborhood of the positive real axis) and has quadratic decay, its derivatives can be 
estimated by 

'K\ 1+p+q p\q\ 



We choose c so large that the following conditions hold, 

ifei 1 

O 16 if, - ~r-<^- (3-22) 

{ujq — o) c 2 

We proceed to prove (]3.20p by induction in I. For I = there is nothing to prove. 
Thus assume that (|3.20p holds for a given I. Using the induction hypothesis together 
with (|3.2ip . we can then estimate the derivatives of the product Wt/j^ as follows, 



Kw>i < e(:)e(i)(t) 



(K\ 1+a+b a! fa! c 1+l+p - a K"- b (p-a)\(q-b)\ 



a=0 v ' b=0 

c i+i+p k 1+ i pi ( , 

a=0 v ~" / 6=0 



A +l+pgl+q p \ q \ P/ K \ " '< 

u 2+l+q l\ Z^\2c 



u 2+b u l+q-b n 
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According to (|3.22p . the two remaining sums can be bounded by the geometric series 
E^=o 2 " m = 2, and thus 



\dPdl(W^)\ < 4 



J.+1+P x l+q pi q\ 

Next we differentiate the integral equation (|3.19l) . 

p 



(3.23) 



r=0 







(v - u) ( e ~ 2l ^ v - u) - 1)1 dP- r (V^ (/) (i>)) dv 



(note that, since Imw < 0, the factor e 2wv gives an exponential decay of the integrand 
as v — > oo). After manipulating the partial derivatives as follows, 



grgq 



e(« - «) ( e - 2 ^(^-«) _ i)l = (_^,)« 



6(« - u) ( 0, ) ( e -«"(«-«) - 1) 



the resulting f-derivatives can all be integrated by parts. The boundary terms drop out, 
and we obtain 



r=0 



-2ioj(v—u) 



i) a 



u — V 



W r (WiP {l) (v)) dv. 



Since u is in the lower half plane, we have the inequality \e 2tul ( v u )| < 1. We conclude 
that 



aja2</.<' +1 >(t0 



£ 2 E 

r=0 



0, 



u — v 



dv . (3.24) 



The ^-derivatives in the curly brackets can act either on one of the factors (u — v) or on 
the function Wip®. Taking into account the combinatorics, we obtain 



r=0 



[v — u) 



ffP-raq-\ 



dv 



Using (|3.23p . we get 



p r 



r=0 s=0 



si (r-s)\l\ 



r—s ( J+l+p-r j^l+q+s 



(v — u) s 



V 



2+l+q+s 



dv . 



Introducing the new variable r = ^, the integral can be computed with iterative integra- 
tions by parts, 



f 

■J u 



1 



(v — u) s 

v 2+l+q+s ' u l+l+q 

l (l + q)\ r 1 



dv 



u i+l+q {i + q + s -)\ J Q 
1 (l + q)\s\ r 1 



[I - TV T l+q dT 



(1 - T Y ^-T l+ 1 +S dT 

y ' dr s 

T l+q+s dT 



1 (l + g)!s! 

u l+l+q (l + l + q + s ) 
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We thus obtain 



< 8 



c i+i+p j£i+q 



V r 

EE 



pi (q + s)l K s r r 



u 1+l +i f^j^ir- s)ll\ (wc) r (l + l + q + s)\ 



Using the elementary estimate 
(q + s)\(l + q)\ 



+ 1 



q + s 



< 



(1 + / + <? + s)! q + l + 1 q + l + 2 q + l + s + l ~ l + l' 



we obtain 

d?d^ l+1 \u) 



< 



,l+l+P K l+ q plq] ^f R ^r r 



u l+l+q a + D) 



The last sum can be estimated by an exponential, 

V T OO 

El / r \r-s ^ 1 f r \ a v-^ 1 / r \ a ( r 

J^syAK) ^^a\\K) ^^al\K) = ^\k 

According to (|3.22p . we can now estimate the remaining sum over r by a geometric series, 

P / Tj r \ r oo 



£(§) exp (^) * £ 



r=0 



j_ \ r 



r=0 



< 2 



We thus obtain 

dldl^ l+l \u) 



c i+i+p K l+q pi q\ c 2 + l +v x q p\ ql 

< 16 ^—r-. — < 



u l+l+q (I + !)! 

where in the last step we again used (13.221) . 



(1 + 1) 



Proof of Theorem\3M According to d3151 13TT8D . 

oo 1 

Expanding ifj^ in a Taylor series in w, we obtain the formal expansion 

oo oo „ 

^(2t(w + C))' ^ 

Lemma 13.41 allows us to estimate this expansion for every to £ Eq n Bg(uJo) with Imw < 
as follows, 




This expansion converges uniformly for |£| < |. Similarly, one can show that the series of 
(^-derivatives also converge uniformly. Hence we can interchange differentiation with sum- 
mation, and a straightforward calculation shows that the Cauchy-Riemann equations are 
satisfied. Thus the above expansion allows us to extend <j> analytically to the ball \Q\ < |. 
Since the constant c is independent of Imw, we thus obtain an analytic extension of </> 
across the real line. ■ 
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3.2 A Continuous Family of Solutions near u = 



In Theorem 13.21 we made no statement about the behavior of the fundamental solutions <f> 
at u = 0. Indeed, we cannot expect the solutions to have a holomorphic extension in 
a neighborhood of u = 0. But at least, after suitable rescaling, these solutions have a 
well-defined limit at uj = 0: 

Theorem 3.5 For every angular momentum number n, there is a real solution of the 
Schrodinger equation {2.1J$ for uj = with the asymptotics 



i 



lim it""* <W«) = with [i := \/A n (0) + ^. (3.25) 

u— >oo ,/7r V 4 



T/iis solution can be obtained as a limit of the solutions from Theorem \3.S\ in the sense 
that for all u € R, 



i 



<fio( u ) = nm ^ 2 <M n ) anc ^ < Poi u ) = nm 2 p( u ) ■ 

Note that the X n are the eigenvalues of the Laplacian on the sphere. They are clearly 
non-negative, and thus the parameter fi in (|3.25p is positive. 

Unfortunately, the function <fo cannot be constructed with the iteration scheme (j3. 16j) 
because if we put in the Green's function for uj = (which is obtained from (|3.14p by 
taking the limit uj —* 0), we get the for </>W the equation 

0«(u) = I (v-u)W(v)dv , 



and since W decays at infinity only quadratically, the integral diverges. To overcome this 
problem, we combine the quadratically decaying part of the potential with the unperturbed 
operator. More precisely, for any uj in the set 

F := {uj G C I Imu < and \uj\ < (16afc) -1 } , 

we write the Schrodinger equation as 

(-^ + ^p-^ 2 )^) = -W(u)<Ku), 

where /j,(uj) = (X n (uj) — 2akuj + |)a . The potential W is continuous in uj and bounded by 

\W(u)\ < for all uj e F . (3.26) 

The solutions of the unperturbed Schrodinger equation can be expressed with Bessel func- 
tions, 

h\{u) = v/™ , h 2 (u) = JT^Ly^uu) . 

They have the following asymptotics, 

h\{u) ~ cos(wn) , h,2(u) ~ sin(u;u) if wm 3> 1 

hl {u) ~ -y^i-^ , w „ r ^ u-^j Jfwu<1 . 

T(/i + 1)2^+5 0FaA4 
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The Green's function can be expressed in terms of the two fundamental solutions by the 
standard formula 

s(u,v) = e(v-u) M«)M«)-M«)M«) 



w(hi,h 2 ) 

where w(hi, h 2 ) = h!^h 2 — h\h' 2 = is the Wronskian. The perturbation series ansatz 

00 

1=1 

now leads to the integral equation 

<j)( l+1 \u) = / S(u,v) W{v)<t> {l \v)dv . (3.28) 



We choose the function <fyP' such that its asymptotics at infinity is a multiple times the 
plane wave e~ tu)U , whereas for u = 0, it has the asymptotics ([3. 25ft . 

^(u) = ^ (hi - ih 2 )(u) . (3.29) 



Lemma 3.6 For any fixed n there is uq S R such that the iteration scheme A3. 291 W. 



converges uniformly for all u > uq and u E F. The functions <p defined by {3.21 ) are 
solutions of the Schrddinger equation {2.1J$ with the asymptotics 



c 

< - 

u 



and a constant c = c(n). 



Proof. Using the asymptotic formulas for the Bessel functions, one sees (similar to the 
estimate [3j eqn. (4.4)] for fj, = I + ^ and integer I) that for all v > u and u 6 F, the 
Green's function is bounded by 



/ „ \ M 1 



| SM | SC ^-(«o^_H_j °(_i_) < 3 ' 30 > 

Similarly, we can bound the Bessel functions in (13.291) to get 

1 < \<f>W\ e~ liauJU ( ^— V < C. (3.31) 

C V 1 + \uj \u J 



Let us show inductively that 



M+l fr<„\ 1 



fCc 



| ^ |so ^-. ( __J (_) . (3.32) 
For Z = there is nothing to prove. The induction step follows from (|3.28l 13.261 13.30j) 

+ J u \l + \w\vj v 3 \ v J 

1+ u uj \ u J L. v 2 
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The lemma now follows immediately from (|3.32l I3.31|) and by differentiating the se- 
ries (|3.27p with respect to u. ■ 

Proof of Theorem \3.5[ ^^From the asymptotics at infinity, it is clear that 

4> if u ^ 
4>q if u! = . 

Denoting the w-dependence of <f> by a subscript, we thus need to prove that for all n£l, 
lim <f} u (u) = cp (u) , lim <jf u (u) = (j>' (u) ■ (3.33) 

To simplify the problem, we first note that for n on a compact intervals, the continuous 
dependence on u follows immediately from the Picard-Lindelof Theorem (i.e. the continu- 
ous dependence of solutions of ODEs on the coefficients and initial values). Thus it suffices 
to prove (13.330 for large n. Furthermore, writing the Schrodinger equation as 

(d u - ito)(d u + iu) § w = -U 4> , 

the potential U has quadratic decay at infinity. Thus, after the substitution {d u — iuS) = 
e lLJU d u e~ luJU , we can multiply the above equation by e~ luJU and integrate to obtain 



oc 



e — «(a u + icu)^(n) = / e~ MV U(v) tf> u (v) dv . 

J u 

Here we emphasized the w-dependence by a subscript; note also that the integral is well- 
defined in view of the asymptotics of (j)^ at infinity. This equation shows that (j)'^ converges 
pointwise once we know that (f>w(u) converges uniformly in n. Hence it remains to show 
that for every e > there is no and 5 > such that for all u> G F with \u>\ < 5, 

\<fru(u) — 4>o( u )\ < £ for all n > no- (3.34) 

To prove (|3.34p we use the uniform convergence of the functions <j>u , (|3.29[) , to choose 5 
such that for all uj £ F with \oj\ < S, 

l#W-0O O) («)l < I for all n> n . 

According to Lemma |3.6| we can by choosing no sufficiently large arrange that 

\^H U ) - 4>uj(u)\ < ^ for all n > no and w G F. 

Now (|3.34p follows immediately from the estimate 

\K ~ M < \<f> u - #| + l# - <PP\ + l4 0) -M- m 
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4 Global Estimates for the Radial Equation 

Let Y\ and Y 2 be two real fundamental solutions of the Schrodinger equation (|2.14p for a 
general real and smooth potential V. Then their Wronskian 

w := Y((u) Y 2 (u) - Y±(u) K,» (4.1) 

is a constant. By flipping the sign of Y 2 , we can always arrange that w < 0. We combine 
the two real solutions into the complex function 

z = Y 1 + iY 2 , 

and denote its polar decomposition by 

z = pe i<fi (4.2) 

with real functions p(u) > and p(u). By linearity, z is a solution of the complex 
Schrodinger equation 

z" = V z . (4.3) 

Note that z has no zeros because at every u at least one of the fundamental solutions does 
not vanish. 

4.1 The Complex Riccati Equation 

We introduce the function y by 

V = j ■ (4-4) 

Since z has no zeros, the function y is smooth. Moreover, it satisfies the complex Riccati 
equation 

y' + y 2 = V . (4.5) 

The fact that the solutions of the complex Riccati equation are smooth will be helpful 
for getting estimates. Conversely, from a solution of the Riccati equation one obtains the 
corresponding solution of the Schrodinger equation by integration, 



log< = y- (4.6) 

J u 

Using (|4.2p in (|4.4p gives separate equations for the amplitude and phase of z, 

p' = p Re y , if' = Im y , 

and integration gives 

logpl^ = [ V Rey (4.7) 

v\l = finy. (4-8) 



u 



Furthermore, the Wronskian (|4.ip gives a simple algebraic relation between p and y. 
Namely, w can be expressed by w = — Im (z z') = p 2 Im y and thus 



w 



P 2 = ■ (4-9) 
im y 

Since p 2 is positive and w is negative, we see that 

Imy(w) > for all u. (4-10) 
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4.2 Invariant Disk Estimates 



We now explain a method for getting estimates for the complex Riccati equation. This 
method was first used in [6] for estimates in the case where the potential is negative 
(Lemma 14. ip . Here we extend the method to the situation when the potential is positive 
(Lemma I4.2|) . For sake of clarity, we develop the method again from the beginning, but 
we point out that the proof of Lemma 14.11 is taken from [6]. Let y(u) be a solution of 
the complex Riccati equation (14. 5|) . We want to estimate the Euclidean distance of y to a 
given curve m(u) = a + i(3 in the complex plane. A direct calculation using (|4.5p gives 

\-^-\y-m\ 2 = (Rey - a) (Rey - a)' + (Imy - (3) (Imy - (3)' 
2 au 

= {Rey -a) [V - (Rey) 2 + (Imy) 2 - a'] - (Imy - f3) [2 Rey Imy + /3'] 

= (Rey -a) [V - (Rey) 2 - (Imy) 2 + 2/3 Imy - a'] + (Rey - a) 2(Imy - (3) Imy 

— (Im y — /3) \J3 + 2a Im y] — (Imy — (3) 2(Rey — a) Imy 
= (Rey - a) [V - (Rey - a) 2 - (Imy - (3) 2 - a 2 + (3 2 - a'] 

-(Imy - f3) [f3' + 2a(3] - 2a ((Rey- a) 2 + (Imy - (3) 2 ) . 

Choosing polar coordinates centered at m, 

y = m + Re ltp , R := \y — m\ , 

we obtain the following differential equation for R, 

R' + 2aR = cos <p [V - R 2 - a 2 + (3 2 - a'] - sin ip [(3' + 2a(3] . (4.11) 

In order to use this equation for estimates, we assume that a is a given function (to 
be determined later). With the abbreviations 

U = V-a 2 -a' and a(u) = exp^2^ a\ (4.12) 

the ODE (|4.1ip can then be written as 

(aR)' = a [U-R 2 + (3 2 ] cos ip - (a /?)' sin ip . 

To further simplify the equation, we want to arrange that the square bracket vanishes. 
If U is negative, this can be achieved by the ansatz 



^vB( T+ i). r= >M(t-^ iu<0) , (4. 13) 

with T > 1 a free function. In the case U > 0, we make similarly the ansatz 



^4(r-i), A=^(T + i) (U>0) <4.,4) 

with a function T > 0. Using (|4.13l B~T4"|) . the ODE (|4.1ip reduces to the simple equation 
(<yR)' = —(a/3)' simp. If we now replace this equation by a strict inequality, 

(aR)' > -(a/3)' simp, (4.15) 
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Imy 




with R a general positive function, the inequality \y — m\ < R will be preserved as u 
increases. In other words, the disk Bn{m) will be an invariant region for the flow of y. In 
the next two lemmas we specify the function T in the cases U < and U > 0, respectively. 
To avoid confusion, we note that it is only a matter of convenience to state the lemmas 
on the interval [0, ii max ]; by translation we can later immediately apply the lemmas on any 
closed interval. 

Lemma 4.1 Let a be a real function on [0, u max \ which is continuous and piecewise C , 
such that the corresponding function U , \J^.12^ , is negative, 

U < on [0, u m J\ . 

For a constant To > 1 we introduce the function T by 

T{u) = T exp Q TF [M log|a 2 [/|) , (4.16) 

define the functions j3 and R by {4-13 ) and set m = a + i(5. If a solution y of the complex 
Riccati equation d^.5|) satisfies at u = the condition 



\y-m\ < R, (4.17) 
then this condition holds for all u £ [0, u max ] (for illustration see Figure^). 
Proof. For e > we set 

T £ (u) = r expQ^"|^|+ e (l-o) (4.18) 

and denote corresponding functions a, R, m, and a by an additional subscript e. Since 
T e (0) = T(0) and lim e \o^e = T, it suffices to show that for all e > the following 
statement holds, 

\y - m e |(0) < R £ (0) =^ \y - m £ \(u) < R £ (u) for all u G [0, n max ] . 
In differential form, we get the sufficient condition 

\y - m £ \{u) = R £ (u) =^ \y - m £ \\u) < R £ (u) . 
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According to (|4.15p . this last condition will be satisfied if 

(a £ R £ )' > \(a £ p £ )'\ . 



(4.19) 



llFiom now on we omit the subscripts e. 

In order to prove (|4.19p , we first use (|4.13l I4.12p to rewrite the functions a (5 and aR 

as 

UV\^U\T+VWU\T^ 



a/3 
aR 

By definition of T £ flUED, 
It follows that 



a 2 U\T- V\aW\T 



-i 



(4.20) 



\a 2 U\' 



+ ee 



{y/WU\T)> 



\a 2 U\ 
-ee~ u {y/\anT\T- 1 ) 



ee- u (^\aHJ\T) 



if \a 2 U\' > 
if \a 2 U\' <0. 



Hence when we differentiate through (|4.20p and set e = 0, either the first or the second 
summand drops out in each equation, and we obtain (aR)' = \<j(3\' ■ If e > 0, an inspection 
of the signs of the additional terms gives (|4.19p , ■ 



Lemma 4.2 Let a be a real function on [0, u mm ] which is continuous and piecewise C , 
such that the corresponding function U, {^.12^ , satisfies on [0, u max ] the conditions 

U > and U' + 4Ua > . (4.21) 

For a constant Tq > we introduce the function T by 



T(n) = T °V^7' (422) 

define the functions (3 and R by ^4-14\ ) and set m = a + i(3. If a solution y of the complex 
Riccati equation j^.5| ) satisfies at u = the condition 

\y — m\ < R , 

then this condition holds for all u G [0, u max ] (see Figure^). Furthermore, 



Rey > a-VU-To^^- . (4.23) 

2a 



i 



Proof. For e > we set 

T £ = T (a 2 U)-t (1-ee- 
Using (|4.141 14. 12 j) we can write the functions aj3 and aR as 



-\ (To" 1 a 2 U (1 - ee-")- 1 - T (1 - ee~ u \ 
1 



2 

aR = - ( '/;, : a 2 Ci\ st "}'••' I < 
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where we again omitted the subscript e. Differentiation gives 

(oR)' > = \ To 1 (a 2 U (1 - ee-")" 1 )' - \ T (l - ee*)' . (4.24) 

According to the second inequality in (|4.2ip . the function o~ 2 U is strictly increasing and 
thus the expression on the right of (|4.24|) is positive for sufficiently small e. Hence (]4. 19|) 
is satisfied. Letting e — > 0, we obtain that the circle Bairn) is invariant. 

In order to prove (|4.23p we note that in the case T < 1 the inequality is obvious because 
even Re y < a — y/U (see Figure |2J) ■ Thus we can assume T > 1, and the estimate 



Rey > a-R > a--^-{T + 2) 

together with (|4.22p gives the claim. ■ 

If the potential V is monotone increasing, by choosing a = we obtain the following 
simple estimate. 

Corollary 4.3 Assume that the potential V is monotone increasing on [0, u max ]. For a 
constant Tq > with Tq > — V(0) we introduce the functions 

/? = H r °-£)' r = K To+ ^)- <425) 



If a solution of the complex Riccati equation |^.5| ) satisfies at u = the condition 
y G {z \z-%P\ < R, Rez,Imz > oj U lz \ z - ^ < y, Rez < o| , 
i/zen ^/i?s condition holds for all u £ [0,u max \ (see Figure^. 

Proof. Choosing a = and (3, T according to f)4.25[) . we know from Lemma 14.11 and 
Lemma 14.21 that the circles \y — m\ < R are invariant. Furthermore, we note that the 
arc A in Figure [3] is the flow line of the equation y' + y 2 = 0, and thus it cannot be crossed 
from the right to the left when V is positive. This gives the result in the case that V has 
no zeros. If V has a zero, the invariant disks in the regions V < and V > coincide at 
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Figure 3: Invariant region estimate for monotone V. 



the zero of V. 



The invariant disk estimates of Lemma 14.11 and Lemma 14.21 can also be used if the 
functions a and U have a discontinuity at some v £ [0,u max ], i.e. 

ot[ : = lim a(u) ^ lira, a (u) =: a r , U\ := \\m.U{u) ^ HmU(u) =: U r . 

u/'v u\v ufv uS^v 

In this case we choose the function T also to be discontinuous at v, 

Ti := lim T(u) ^ lim T(u) =: T r , 

u /x; u\v 

in such a way that the circle corresponding to (a r ,U r ,T r ) contains that corresponding 
to (ai,Ui,Ti) (see Figure HJ). In the next lemma we give sufficient "jump conditions" for 
this "matching." 

Lemma 4.4 (matching of invariant disks) Suppose that U\ < 0. Depending on the 
sign of U r , we set 



Tr 



Ti 



(a r -ai) 2 + \Ui + U r 



if U r < 



VWlWr 



(4.26) 
(4.27) 



Let By r be the disks with centers my r = ay r + ij3y r and radii Ry r as given by ^4.13 ) 
or (Jlfl) . Then B t C B r . 

Proof. We must satisfy the condition R r > \m r — mi\ + Ri. Taking squares, we obtain 
the equivalent conditions R r > Ri and 

{Rr-Ri? > (a r - ai ) 2 + ((3 r - ft) 2 . 

This last condition can also be written as 

(a r - ai ) 2 + ((3? - Rl) + (P 2 - R 2 ) < 2 {PiPr - R h Rr) . (4.28) 
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Figure 4: Matching of invariant disks in the cases U r < (left) and U r > (right). 



In the case U r < 0, we can substitute the ansatz (|4.13p into (|4.28p to obtain the 
equivalent inequality 

{a r - ai ) 2 + \U l \ + \U r \ < y/piWUr] + ■ 

Dropping the last summand on the right and solving for T r , we obtain (|4.26p . which is 
thus a sufficient condition. 

In the case U r > 0, we substitute (j4.13|. 14. 14f) into (14.280 to obtain the equivalent 
condition 

(Or ~ atf + \U t \ - U r < VWWr ~ ^) ■ 

Using the inequality \U[\ — U r < \Ui + U r \, replacing the factor Ti/T r on the right by one 
and solving for T r , we obtain the sufficient condition (|4.27p . ■ 



4.3 Bounds for the Wronskian and the Fundamental Solutions 

We now consider the solutions <j> and 4> as defined in Section 13.11 for oj on the real axis and 
set 

. i' . i>' 

»-■?• y -~y 

We keep k fixed. Since taking the complex conjugate of the separated wave equation flips 
the sign of k, we may assume that k > 0. Then uq as defined by (I2.12p is negative. 

Proposition 4.5 If to [u;o,0] ; the Wronskian w(<j),(j)) is non-zero. 

Proof. According to (|2. 13|) . co and 0, have the same sign. From (|4.10p we know that 
the functions y and y both stay either in the upper or lower half plane. In view of the 
asymptotics (|3.1l 13 ■ 2|) . we know that they must be in be in opposite half planes. Thus 

w(J),4>) = (f>4>(y -y) ^ . * 

In the case u £ (u>o,0), we need the following global estimate for large A. 
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Proposition 4.6 For any u\ G R there are constants c, Xq > suc/i that 

Q 

< — for all A > Ao, u G (wo,0), u < Hi . 
A 

The remainder of this section is devoted to the proof of this proposition. Let u\ G R 
and lo G (wo,0). Possibly by increasing iti and Ao we can clearly arrange that V is 
monotone decreasing on [u±,oo). Then we have the following estimate. 

Lemma 4.7 The functions <j> and y satisfy the inequalities 

> 1, Rey(u) < \co\ on[u\, oo). 

Proof. From the asymptotics (|3.2p we know that lim^^oo y{u) = —iuj. Thus for v 
sufficiently large, \y(v) — i\u>\\ < e, and we can apply Corollary 14.31 on the interval [u±,v] 
backwards in u with To = \lo\ + 2e. Since e can be chosen arbitrarily small, we conclude 
that Corollary 14.31 applies even on [1x1,00) with To = \uj\. This means that 

< Imy < |ct»| , Rey < \u>\ on [1x1,00) . 

Finally, we use (|4,9p with w = i\u>\. ■ 




We now come to the estimates for <fi, which are more difficult because we need a 
stronger result. The next lemma specifies the behavior of the potential on (— oo,2ui]. 

Lemma 4.8 For any u\ £ R there are constants c, Ao such that the potential V has for 
all co G (ojq, 0) and all X > Xq the following properties. There are unique points U- < uq < 
u + < u\ such that 

V(u-) = ~, V(u ) = 0, V(u + ) = Q 2 . 

V is monotone increasing on (—00,114.]. Furthermore, 

< c (4.29) 
> log n 2 - log A - c (4.30) 

< ^\V\ on [u+ ,2m], (4.31) 
with 7 as in Ii3.9\ ). 

Proof. We expand V in a Taylor series around the event horizon, 

V = -n 2 + (X + c )(r-r 1 )+XO((r-r 1 ) 2 ). 

Hence for sufficiently large Ao there are near the event horizon unique points u_, no, u+ 
where the potential has the required value. Integrating (12. ip we get near the event horizon 
the asymptotic formula 

u ~ — log(r — n) . 



U + — U- 

711+ 

V'\3 + \ V 2 



22 



Getting asymptotic expansions for u± we immediately obtain (|4.29l FOUj) . Furthermore, 
using (|2.ip to transform r-derivatives into u-derivatives, we obtain in the region (ri,ri + 
e) H (u + , oo) the estimates 

-e~< u < V(u) < Ace 7 " 
c 

\V'(u)\ + \V"(u)\ < Ace 7 ", 

uniformly in A and oj. Hence for sufficiently large Ao, (|4.3ip will be satisified near the 
event horizon. 

In the region r > r\ + e away from the event horizon, V is strictly positive, V > A/c, 
and since the derivatives of V can clearly be bounded by \V'\ + \V"\ < cA, it follows 
that (|4.3ip is again satisfied. ■ 



First we apply Corollary 14.31 on the interval (— oo,u_) to obtain the following result. 
Corollary 4.9 There is a constant c > such that for all uj £ (uq,0) and A > Ao, 

— < Imy < U , |Rey| < — on (— oo, u_\ . 

Also, at u = u- we have an invariant disk with 

ai = 0, Ui = -— , Ti = \/2 . (4.32) 

On the interval we use the method described in the next lemma. 

Lemma 4.10 Assume that the potential V is monotone increasing on [0, u max \. We set 



a = v / max(2 V(u max ),0) 

and introduce for a given constant Tq > 1 the functions U, a, f3, R, and T by \4-lS\ \4-14\ ) 
and 

T(u) = T e 2 - . (4.33) 

If a solution y of the complex Riccati equation |^.5| ) satisfies at u = the condition 

\y — m\ < R , 
then this condition holds for all u £ [0, u max ] . 

Proof. By definition of a, the function U = V — a 2 is negative and monotone increasing. 
Using furthermore that a = e 2au , we can estimate the total variation in (|4.16p as follows, 

TV [M log \a 2 U\ = £(^4a-^ = 4au + log \U(0)\ - log \U(u)\. 

This gives ItPB]) . ■ 

Thus we match the invariant disk ()4.32p to a disk with U r = V(uJ) — a 2 and a r = y2Q. 
From (|4.29p we see that (n + — v,-) a is uniformly bounded, and thus we obtain the following 
estimate. 
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Corollary 4.11 There is a constant c > such that for all co £ (^OjO) and A > Xq, 

— < Imy < cQ. , |Rey| < cQ. on[u-,u+]. 
c 

At u = u + we get an invariant disk with 

< ai < ft, -cUi = -n 2 , Ti < c. (4.34) 

In the remaining interval [n + ,2tii] an approximate solution of the Schrodinger equa- 
tion (|2.14p is available from semi-classical analysis: the WKB wave function 

(j)(u) = V-\ exp (^j VV^j . 

The corresponding function y is given by 

y(u) = W--. 

In order to get an invariant disk estimate which quantifies the exponential increase of (p, 
we choose a such that it also becomes large as V ^> 0. For technical simplicity, we choose 

a (u) = 7 -VW), (4.35) 
giving rise to the following general result. 

Lemma 4.12 Assume that the potential V is positive on [0,u max ] and that 

\V'(u)\ < \ v{u)l, \V"(u)\ < \ v{uf. (4.36) 



We introduce for a given constant Tq > the functions a, U , a, (3, R, and T by {1^.35 



\4-l%\ \4-M\ \4-%ty - If a solution y of the complex Riccati equation (f^.5| ) satisfies at u = 
the condition 

\y — m\ < R , 

then this condition holds for all u £ [0, u, nax ] . Furthermore, 

Rey > - ^\n\. (4.37) 

Proof. A short calculation yields 



64 16 y/v 

TTl . TT 105 3 83 TA/ 7 V' 2 7 V' 
U' + AaU = V2 -V + 



128 64 32 V i 16 W ' 

Using (|4.36p we obtain the estimates 



v v yl 

< u < and U' + AaU > — 

64 2 16 
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Hence the conditions (|4.21|) are satisified, and Lemma 14.21 applies. The inequality (|4.37p 
follows from (|4.23p . the just-derived upper bound for U and the fact that a > 1. ■ 



Matching the invariant disk (|4.34p to the invariant disk with a r = a{u r ) and U r 
V(u r ) — a 2 (u r ) — a'{u r ) with a according to (|4.35p . we obtain 



U r < Q 2 , T r < c. 
We can then apply the last lemma on the interval [u+, 2u\]. 



(4.38) 



Proof of Proposition \4~1>[ Suppose that u < u\. Using the definition of y and y, we can 
rewrite the Wronskian as 

w(4>,4>) = <j></> (y-y) . 
Applying Lemma 14.71 at u = 2u\ gives 



< 



0(2«i) 



1 



Re y(2ui) — \u)\ 



We combine (|Q71) with T = T r satisfying (g^Sj) to get 

Re y > — cCl . 



(4.39) 



(4.40) 



Since the potential V is strictly positive on the interval [tti,2tii], we can, possibly by 
increasing Aq and c, arrange that 



and thus also that 



w > 



Rey > 



16 c 



on [m, 2ui] 



on [u\, 2ui] . 



This inequality allows us to bound the fraction in (g7 



<f>(u) 



W(<j>,(j)) 



< 



cp{2 Ul ) 



(4.41) 
(4.42) 

(4.43) 



Thus it remains to control the last quotient. We omit the accent and use the nota- 
tion p = \<j)\. In the case u < «+, we can use (|4.9p . 



p{u) 



1my(u+) 



/3(u+) 2 Imy(u) ' 

and the last quotient is controlled from above and below by Corollary 14.91 and Corol- 
lary 14.111 Hence, rewriting the quotient on the right of (I4.43P as 



p(u) 



p(u) p(u + ) 



p(2m) p(u + ) p{2u{) 
it remains to consider the case u > u + . Applying (|4.7p and ()4.40p . we obtain 



A := loe 



<P(2 Ul ) 



2«i 



Re y(u) < cO, (2ui — u + ) 



1 



2ui 
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Now we use (|4.3U|) and the fact that the function log is bounded, 



A < c log A / W 



2u 



Estimating the last summand with (j4.41H . 



2u\ r2ui /T 

y/V > I V7 > 



we conclude that for large A this summand dominates the term clog A, and thus (|4.43|) 
decays in A even like exp(— y/X/c). ■ 



5 Contour Deformations to the Real Axis 

In this section we fix the angular momentum number k throughout and omit the angular 
variable ip. We can again assume without loss of generality that k > 0. Also, since 
here we are interested in the situation only locally in u, we evaluate weakly. Thus we 
write the integral representation (|2.16p for compactly supported initial data ^>q an d a test 





function ^^(Ix S 2 ) 2 as 



<ri, tt(t)> = o~~~- Y] I™ ( I ~ I ) d " e~ iuJt <r?,Qn(w) SooM *o> • (5.1) 

The integration contour in (|5,ip can be moved to the real axis provided that the integrand is 
continuous. In the next lemma we specify when this is the case and simplify the integrand. 
For lo real, the complex conjugates of and <fi are again solutions of the ODE. Thus, apart 
from the exceptional cases lo G {0,ujo}, we can express </> as a linear combination of <j) 

and cj), 

= a(j> + p (j) (w£i \ {0, ll>o}) . (5.2) 

The complex coefficients a and j3 are called transmission coefficients. The Wronskian of cf> 
and 4> can then be expressed by 

w(j>,4>) = f3w(4>,<f>) = 2inp, (5.3) 

where in the last step we used the asymptotics ()3. 1 j) . Furthermore, it is convenient to 
introduce the real fundamental solutions 

(f>i = Re cj) , 4>2 = Im <j> , 

and to denote the corresponding solutions of the wave equation in Hamiltonian form 
by 

Lemma 5.1 // the Wronskian w(<j>, <j)) is non-zero at lo € M \ {0,loq}, then the integrand 
in H5.1\) is continuous at lo and 



(lhn-lun)(Q n (u; + ie)S oo (L0 + i e)^)(r^) = — - £ *r , (5.4) 

e/0 e\0 LOU 

a, 6=1 
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where the coefficients t a b are given by 

in = 1 + Re-, ti2 = t 2 i = _Im ^> *22 = 1-Re— . (5.5) 

Proof. We start from the explicit formula for the operator product Q n Soo given in [5l 
Proposition 5.4]. Since the angular operator Q n (u>+ie) can be diagonalized for e sufficiently 
small, the kernel g(u, u') is simply the Green's function of the radial ODE, i.e. for cu in 
the lower half plane, 

9{u > u) - W (U) \ku)kw) ifuw. 1 (5 - 6) 

whereas the formula in the upper half plane is obtained by complex conjugation. Using 
that lim e yo 4> = nm e\o 4> an d lim e yo 4> = nm e\o k we find that 

( lim — lim ) g(u, u) = 2i Im g(u, u) , 

e/0 e\0 



and a short calculation using (|5.2l 15 .3|) gives 

2 



(lim- lim = -1 ^ U0» /(u') 



a,6=l 

with t a fe according to (|5.5p . 

Except for the function g(u,u'), all the functions appearing in the formula for Q n S 
in Proposition 5.4] are continuous on the real axis. A direct calculation shows that 

( lim - lim )(Qk, n ( UJ + i £ ) £oo(w + is) *&) 

e/O e\0 



.nE^ 1 o ij*>^(*.) 

a,6=l v 7 



with c?/u given by (|2.7|) . Since the Vfj, are eigenfunctions of the Hamiltonian, we know 
according to (|2.3|) that A^^ = (u — Using furthermore that the operator A is 

symmetric on L 2 (dn), we conclude that 

where in the last step we used (|2.8j) . ■ 



Let us now consider for which values of w and n the contour can be moved to the real 
axis. According to Proposition 14.51 the Wronskian w((f),(j)) is non-zero unless to £ [ujq,0]. 
We now analyze carefully the exceptional cases lo = 0, loq. From Theorem l3-H Theorem l3.2l 
and Theorem [33] we know that the functions 4> and <f> w = uj^cj) are continuous for all u € R. 
If ak = and u = 0, the functions <j> and ^ degenerate to real solutions with the 
asymptotics 

lim <p(u) = 1 , lim u^~z <fo(tt) = \— ■ 
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Figure 5: The integration contour D$- 



Noting that the function 

is monotone increasing, the potential V, (|2.15p . is everywhere positive. Hence solutions 
of the Schrodinger equation (|2.14p are convex. This implies that the functions 4> and 
do not coincide, and thus their Wronskian is non-zero. As a consequence, the Green's 
function (|5.6p . and thus the whole integrand in (|5.ip . is bounded and continuous near u = 
(note that (15. 6p is invariant under rescalings of 0, and thus we can in this formula replace <j> 
by <t>w)- in the case ak ^ and uj = 0, the function 0q is real, whereas <fi is complex, and 
thus w(4>, 4>o) 7^ 0. If on the other hand ak ^ and a; = ujq, (j> is real and <p is complex, 
and again w(<p,cpo) / 0. Hence the integrand in (|5.ip is continuous and bounded at the 
points u; = 0,wo- We conclude that for every n £ N, the integrand in (|5,ip is continuous 
on an open neighborhood of on uj £ R\ (cjq, 0). Furthermore, according to Proposition 14.61 
w((j),(j)) ^ if uj £ (ujo,0) and A is sufficiently large. We have thus proved the following 
result. 

Proposition 5.2 There is 5 > and uq G N suc/i that for every £ Cg°(M x S 2 ) 2 , the 
completeness relation 

+ E f ^ n 500 dw 

n<n ^ D <5 

holds, with the contour Dg as in Figure 0. 

We point out that the contour D$ passes along the line segment [ujq,ujq + 8) twice, once 
as the limit of the contour in the lower half plane, and once as limit of the contour in the 
upper half plane. These two integrals can be combined to one integral over [ujq,ujq + 5) 
with the integrand given by (15.41) . 

Let us now consider how the remaining contour integrals over C £ can be moved to the 
real line. According to Theorems 13.11 and 13.21 the functions cj) and <p have for every n < uq 
and for every uj £ (uj q,0) a holomorphic extension to a neighborhood of uj. Thus their 
Wronskian is also holomorphic in this neighborhood, and consequently they can have only 
isolated zeros of finite order. Since w(4>,4>) ^ for ui near and u>o, we conclude that the 
numbers of zeros must be finite. Since we only need to consider a finite number of angular 
momentum modes, there is at most a finite number of points uj%, . . . , ujk £ (^o, 0), K > 0, 
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where any of the Wronskians w(4> n , <j) n ) has a zero. We denote the maximum of the orders 
of these zeros at U{ by k £ N. 

The above zeros of the Wronskian lead to poles in the integrand of (|5.ip and correspond 
to radiant modes. We will prove in Section [7] by contradiction that these radiant modes 
are actually absent. Therefore, we now make the assumption that there are radiant modes, 
i.e. that the Wronskians w(<fi n , <p n ) have at least one zero on the real axis. As a preparation 
for the analysis of Section [7J we now choose a special configuration where radiant modes 
appear, but in the simplest possible way. We choose new initial data 

$o = V{H) , (5.7) 

where V is the polynomial 

K 

V(x) = lo(oj-ujq) {x-Uif 1 ' 1 Y\[x-uJi) k . 

i=2 

Then $o again has compact support, and using the spectral calculus, the corresponding 
solution $(i) of the Cauchy problem is obtained from (|5.ip by multiplying the integrand 
by V(u>). Then the poles of the integrand at u>2, ■ ■ ■ ,wk disappear, and at w\ a simple 
pole remains. Subtracting this pole, the integrand becomes analytic, whereas for the pole 
itself we get a contour integral which can be computed with residues. 

Let us summarize the result of the above construction with a compact notation. For 
a test function rj £ C^°(M x S 2 ) 2 we introduce the vectors •rf }n by 



Proposition 5.3 Assume that there are radiant modes, K > 1. Then the Cauchy devel- 
opment <3?(t) of the initial data |5. 7| ) satisfies the relation 



1 f°° 

<r],$(t)> = — V / e~ iujt (■n ujn ,T ujn <5> ujn ) C 2duj + e~ iuJlt V (rf 1 "', a n ) C 2 
neN 00 n<n 



Here uj% £ (uo,0). The (o" n ) n= i v .. ino are vectors in C 2 , at least one of which is non-zero. 
The matrices T^ n have the following properties, 

(1) Ifujg [uj ,0] orn> n , 

(T" n U = C (5-8) 



with t a fe according to l[5.5\ ). 
(2) For each n, the function T wn is continuous in to £ M. and analytic in (wo,0). 



6 Energy Splitting Estimates 

In this section, we consider the family of test functions 

t)l(u) = T]{u + L) 

for a fixed 77 £ Cq°(1R x S 2 ) 2 . Our goal is to control the inner product <?7l, &(t)> m the 
limit L — > 00 when the support of r/x moves towards the event horizon. Our method is to 
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split up the inner product into a positive and an indefinite part. Once the indefinite part 
is bounded using the ODE estimates of Section [H we can use the Schwarz inequality and 
energy conservation to also control the positive part. 

We choose u\ G R and general test functions ?),( 6 Co°((— oo,ui) x S 2 ) 2 which are 
supported to the left of u% (for later use we often work more generally with £ instead 
of ^o)- Since for each fixed n, the T wn are continuous and the eigensolutions ^>^ n (u) are, 
according to Theorem 13 .1\ also continuous in u, uniformly for u £ (—00, u±), we have no 
difficulty controlling the expressions {rf n , T um C l 0}n )^ for n < uq and uj E [wo,0]. Hence 
we only need to consider the case when the matrix T wn is given by f)5.8[) . Using (|5.5p . the 
eigenvalues X± of this matrix are 



1 ± 



(6.1) 



In order to determine the sign of these eigenvalues, we first use the asymptotics (|3.1l 13.21) 

v \ / ~~7 

to compute the Wronskians w(4>,<p) = —2iuj and w((j>,(p) = 2iVL. Furthermore, we obtain 
from (|5.2p and its complex conjugate that 

w(U) = (\a\ 2 - \P\ 2 ) w(<pj) . 
Combining these identities, we find that 

12 iff - 



uJ 

o 



(6.2) 



From (|6.1l I6.2|) we see that in the case uj [cjq,0], where u and £1 have the same sign, 
the eigenvalues X± are both positive. However, if u G (uq,0), one of the eigenvalues is 
negative. This result is not surprising, because the lack of positivity corresponds to the 
fact that for uj € [cc>o,0] the energy density can be negative inside the ergosphere. In the 
case when T^ n is not positive, we decompose it into the difference of two positive matrices, 



for uj G (loq, 0), n > no , 



where 



In the next lemma we bound the integral over Tf 1 using ODE techniques. 
Lemma 6.1 For any e > we can, possibly by increasing uq, arrange that for all L < 0, 

rO 



[ \{riT^ n C n )cA du, < e. 



n>n J ^0 

Proof. Using f\6.2\ 15 .3|) we can estimate the norm of T_ by 



lA. 



Iff 



< 



1 



2\uQ\ 



Hence 



n>n 



IC 2 



duo < 2 V 



n>ng 



2\ff 



\w{4>,4>)\ 2 



u \w((j),(j))\ \w(<j>,< 



Ljill dui . 
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Writing out the energy scalar product using [SJ eq. (2.14)] and expressing the fundamental 
solutions VE^ 71 in terms of the radial solution (ft, one sees that 

\nt n \ < csu P \ VL 4>\ , \cn < csup|c<AI , 

R K 

where the constant c = c(uj) is independent of A. Now we apply Proposition 14,61 and use 
that the eigenvalues A n grow quadratically in n, (|2.1ip . ■ 



Lemma 6.2 There is a constant C > such that for all L > 0, 

E / \(vz n ,T" n cnc*\<^ < c. 

n6N MN,o] 

Proof. First of all, using the positivity of the matrix T, 

neN Mk>,o] 

< ^E/ ((crc)c^(rrnc^. 

^ neN MN,o] 

The two summands can be treated in exactly the same way; we treat the summand 
involving rffj 1 because of the additional L-dependence. Applying Proposition 15.21 and 
dropping all negative terms, we get 



/ {r,f'\ T^vDc* dco < < VLl H (H - u ) i lL > 

Jm\[u} ,o] 

+ E f (VL n ,T? n VL n )c*<^ + E / l <7 ^> (QnS^HiH-^VL)^ 



Using the asymptotic form of the energy scalar product and the Hamiltonian near the 
event horizon, it is obvious that the first term stays bounded as L — * oo. The second 
term is bounded according to Lemma 16.11 For the contour integrals we can use the for- 
mula (15. 4p on the real interval [ujq,ujq + 5). Since Theorem 13.11 gives us control of the 
asymptotics of fundamental solution <j) uniformly as u — > — oo, it is clear that the integral 
over [ujo,ujq + 5) is bounded uniformly in L. For the contour in the complex plane, we 
cannot work with (|5,4p . but we must instead consider the formula for the operator prod- 
uct Q n Soo give in [5J Proposition 5.4] together with the estimate for the Green's function 
given in Lemma 16.31 below. ■ 



Lemma 6.3 For every u £ D$ with uo ^ u>q, there are constants C, e > and tioGR such 
the Green's function satisfies for all u E D$ n B t (ui) the inequality 



\g(u,v)\ < C for all u,v < uq . 
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Proof. It suffices to consider the case Imw < 0, because the Green's function in the upper 
half plane is obtained simply by complex conjugation. By symmetry, we can furthermore 
assume that u < v. Thus, according to (|5.6p . we must prove the inequality 



4>{v) 



W(<j>,(j)) 



< C 



for all u < v < uq 



According to Whiting's mode stability [TT], the Wronskian w(cf>, </>) has no zeros away from 
the real line, and thus by choosing S so small that B${uj) lies entirely in the lower half 
plane, we can arrange that \w(4>, 4>)\ is bounded away from zero on B${Co). Hence our task 
is to bound the factor \<fi(u) 4>(v)\. Solving the defining equation for w(4>,(p) for cj) and 
integrating, we obtain 



-w(4>,<j)) 



Substituting the identity 
1 



-2iQu 



d 



(j){u) 



2iQ. du 



~,2iQu 



4>(uf 



du 

4>(u)' 



1 d 
2iQ du 



the integral over the last term gives a boundary term, 



" 1 d 
2iQ du 



1 



1 



<j)(uy J 2in 0( u )2 

The integral over the other term can be estimated by 



(f){uf 



1 

2tt 



PU 






e -2iflu 

















dv < 



-2Im O v 



n 



(e- inu 4>(u)Y 



dv 



Using the asymptotics (|3.ip one sees that the last integrand vanishes at the event horizon. 
From the series expansion for <fi, <\3.7\ I3.8p . we see that this integrand decays even expo- 
nentially fast. Therefore, the last integral is finite, uniformly in v and locally uniformly 
in to. Collecting all the obtained terms and using the known asymptotics (|3.ip of 4>, the 
result follows. ■ 



Lemma 6.4 For any e > we can, possibly by increasing uq, arrange that for all L > 0, 



e . 



n>no ' 

Proof. Again using positivity, it suffices to bound the terms 

£ fivT^rvDc^du; and £ /V", OTc^ • 

n>n n>no Ju) 
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They can be treated similarly, consider for example the first term. For any n\ > uq, 

< <At}l, H(H- cjo) A Vl > + X) / ((%)™,r(%)™) c2 dw 
+ X) f l<-4^ (<9n -Soo H(H- u) Q ) Arj L )> | dw . 

n<n D{ 

Here .A is the angular operator. When it acts on a test function, we always get rid of 
the time-derivatives with the replacement idt — > H. We now argue as in the proof of 
Lemma 16.21 (with ijl replaced by Atjl) and choose n\ sufficiently large. ■ 



7 An Integral Representation on the Real Axis 

We now use a causality argument together with the estimates of the previous section to 
show that the radiant modes in Proposition 15.31 must be absent. This will be a contradic- 
tion to the assumption that there are radiant modes, ruling out the possibility that there 
are radiant modes at all. This will lead us to an integral representation of the propagator 
on the real axis. 

Let us return to the setting of Proposition 15.31 Choosing the ^-dependence of r] 
such that it is orthogonal to the angular wave functions (^a in )n<n except for one n, 
and choosing the u-dependence of r\ such that it is orthogonal only to one of the plane 
waves gi^i-^o)^ we can dearly arrange that 

limsup|a(L)| =: k > where a(L) := ^ <if L in ,a n > C 2 . (7.1) 

L ^°° n<n 

Furthermore, we choose t]l such that its support lies to the left of supp<3?o, i-e. 

dist(suppr/^, supp $0) > L for all L > 0. 
Due to the finite propagation speed (which in the (i, u)-coordinates is equal to one), 

suppr/L n supp$(t) = if |*| < L. 

Hence for all L > 0, 

2L 



<» - 7TT / e iuJlt <VL, *(*)> 



1 f°° sm (( u ~ Ui)L) rr^n^wnx j. , „( T \ 



neN' 

We apply Lemma 16. II and Lemma 16.41 with e = k/(8it) to obtain 

»o 

Kn-r'.'l'-' 11 ^") -A ,L: < 



n>n Ju ° 
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Furthermore, Lemma 16.21 gives rise to the estimate 



f sin((u; — loi)L) 
W\[u o ,0] (w-Wi)L 



L ! 



IC- 



du 



< C sup 

K\N,0] 



sin((u; — u\)L) 



(u — oj{)L 



and since this supremum tends to zero as L — > oo, we conclude that the expression on the 
left vanishes in the limit L —* oo. Combining these estimates with (I7.ip . we obtain 



lim sup 



E 



^in^-^L) {rit n^n^n )c2duJ 



n<n JuJ 



(oj — OJ\)L 



> TTK . 



(7.2) 



Since the matrices T un are continuous in u and the fundamental solutions ^>^ n (u) are 
according to Theorem 13.11 uniformly bounded as u — * — oo, there is a constant C such that 

|(C\T wn $nc2| < C for all L > and w G (u; ,0), n < n . 

Hence we can apply Lebesgue's dominated convergence theorem on the left of (|7.2|) and 
take the limit L — > oo inside the integral, giving zero. This is a contradiction. 

Since radiant modes have been ruled out, we know that the Wronskian w(4>, 0) has no 
zeros on the real axis. Thus we can move all contours up to the real axis. This gives the 
following integral representation for the propagator. 

Theorem 7.1 For any initial data £ C(J°(lRx S 2 ) 2 , the solution of the Cauchy problem 
has the integral representation 



fceZ nelN- 7 00 a,b=l 



wft t/ie coefficients t a b as given by A5.51 \5.2\i . Here the sums and the integrals converge 



8 Proof of Decay 

As was recently pointed out to us by Thierry Daude, there is an error on the last page of the 
printed version of this paper (which is identical with the previous online version). In order 
to keep the best possible agreement between the printed and online versions, in this section 
we simply give the erratum which corrects the printed version. The error was that the 
inequality (8.3) in the printed version of this paper cannot be applied to the function ^f(t) 
because it does not satisfy the correct boundary conditions. This invalidates the last two 
inequalities of the printed version, and thus the proof of decay is incomplete. We here fill 
the gap using a different method. At the same time, we will clarify in which sense the 
sum over the angular momentum modes converges in [5J Theorem 1.1] and Theorem \7.1\ 
an issue which so far was not treated in sufficient detail. The arguments so far certainly 
yield weak convergence in L? ; here we will prove strong convergence. 

Our method here is to split the wave function into the high and low energy components. 
For the high energy component, we show that the L 2 -norm of the wave function can be 
bounded by the energy integral, even though the energy density need not be everywhere 
positive (Section 18. lh . For the low energy component we refine our ODE techniques 
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(Section 18. 2p . Combining these arguments with a Sobolev estimate and the Riemann- 
Lebesgue lemma completes the proof (Section 18. 3[) . 

We begin by considering the integral representation of Theorem 17.11 for fixed k and a 
finite number no of angular momentum modes, 

poo 

$ n °(i,r,tf) := / duj e~ iu)t § n °{uj,r,$) , (8.1) 



where for notational convenience we have omitted the (/^-dependence (i.e. the factor e lk ^ : 
and <& n ° (uj) is defined by 

, 1 no 2 

n=l a,6=l 

(as in [5], we always denote the scalar wave function by <£, whereas \& = (<£, c^) is a two- 
component vector). We recall that for large uj, the WKB-estimates of [5j Section 6] ensure 
that the fundamental solutions &\ un go over to plane waves, and thus, since the initial 
data is smooth and compactly supported, the function $ n °(o;, r, •&) decays rapidly in uj 
(for details on this method see [8l proof of Theorem 6.5]). As a consequence, & n ° and its 
derivatives are, for r and $ in any compact set, uniformly bounded in time. Our goal is 
to obtain estimates uniform in no, justifying that, as no — > oo, <3?"° converges in L 2 oc to 
the solution of the wave equation. 

To arrange the energy splitting we choose for a given parameter J > 0, a positive 
smooth function Xh+ supported on (J, oo) with Xh+|[2J,oo) = 1- We define Xh- by Xh-(^) = 
Xh+{— and set Xl = 1 — Xh+ — Xh-- We introduce the high-energy contributions <3E>2± 
and the low-energy contribution <E>l° by 

Ki(uj,r,ti) = X H±H<l no (^r-,tf), $^{u,r,6) = Xl M l> n ° (w, r, 0) . 
8.1 L 2 -Estimates of the High-Energy Contribution 

We recall from [5j (2.5)] that the energy density of a wave function $ in the Kerr geometry 
is given by 

£($) = ( ^ + ° 2 ) 2 - a 2 gin 2 ^ 1^,2 + A ^2 

+ sin 2 i? |9 cos ^| 2 + (— ^- - ^ A: 2 |$| 2 . (8.2) 
\ sin v Ay 

Note that the energy density need not be positive due to the last term. However, the next 
theorem shows that the energy integral 



/oo rl 
dr J dcostf £($(*)) 



(which is independent of time due to energy conservation), in the high-energy region is 
both positive and can be bounded from below by the L 2 -norm. In what follows, we only 

-Hi- 



consider $2+ because $2- can be treated similarly. 
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Theorem 8.1 There exists a positive constant Jq (depending only on k, but independent 
of uq and *$>o), such that for all J > Jq the following inequality holds for every t: 



j2 poo ( r 2 , 2\2 pi 

Em + ) > y I a dr I ^«»* wc*)i a - 

The remainder of this section is devoted to the proof of Theorem 18.11 We begin with 
the following lemma. 

Lemma 8.2 Let 5,<3? be measurable functions, g real and <I> complex, such that <I> and g<& 
are in L 1 (R). Then 



du duj' mm(g(uj),g(uj')\ $(a/) > inf<7 

Jr Jr ^ ' 



2 



Proof. Using a standard approximation argument, it suffices to consider the case that g 
and $ are simple functions of the form 

A A 

= ^9a X{K a ) , $M =^2^a X(K a ) , 

a=l a=l 

where x{K a ) is the characteristic function of the set K a , and (Kg)a=i,...,A forms a partition 
of R. Then the above inequality reduces to 

A A 

mm(g a ,g b ) $ a \K a \W b \K b \ > ming ^ <f> a \K a \W b \K b \ . 
a,b=l a, 6=1 

In the case A = 2 and g\ < g2, this inequality follows immediately from the calculation 

9i | ci | 2 + 5i (ci c^+clc 2 ) + 52 |c 2 | 2 > gi \ci\ 2 +gi (ci c^ + 'cjc 2 ) + g\ |c 2 | 2 = gi\ci + c 2 \ 2 , 

where c a := $ a |if a |- in the case A = 3 and 51 < 52 < 53, we get 

51 (|ci| 2 + 2Re(ciC2" + ciq)) + 52 (|c 2 | 2 + 2Re(c 2 Q)) + 53 |c 3 | 2 

> 5i (|ci| 2 + 2Re(ciC2" + ciq)) + 52 |c 2 + c 3 | 2 > 51 |ci + c 2 + c 3 | 2 . 

The general case is similar. ■ 



The next lemma bounds the L 2 -norm of <£h+ an d its partial derivatives by a constant 
depending on no and \&o- 

Lemma 8.3 There is a constant C = C(uq, ^q) such that for every t, 



2 1 „2\2 /■! , 3 



(r 2 + a 2 ) 



dr f dcos^(|^(t)| 2 + |^Ui)| 2 + El^ + W| 2 ) <C 
J ~ x k=i 



Proof. It suffices to consider one angular momentum mode. For notational simplicity we 
omit the angular dependence. Since 3>2+ an d its derivatives are locally pointwise bounded 
uniformly in time, it follows that their L 2 -norms on any compact set are bounded in time. 
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Near the event horizon, we work with the fundamental solution (ft in the Regge- Wheeler 
variable u (see (2.18, 5.2)] and Section 13.11 Then for any sufficiently small Uq, the 
integral of \&\ 2 over the region u < uq can be written as 



u 

du\(j)\ 2 



where it is now convenient to write our integral representation (|8.1|) in the form 
4>{t,u) = I dw (h+(uj) <f> w (u) + h-(u) <j*u(u)^ 



Here the functions h± have rapid decay and, as they are supported away from the 
set {0,ujo}, they are also smooth (see Section IXTj) . Using the Jost representation (|3.71 
3.8j> . the function (f>(t,u) can be decomposed as 



4>(t,u) = <fr+(t,u) + 4>-{t,u) + p(t,u) 

where 

/oo 
-oo 

\p(t, u)\ < C e 7 " for all t, where C, 7 > 0. 

Note that the smoothness of h± implies that <j>± decay rapidly in u. From the exponential 
decay of the factor e yu it is obvious that the L 2 -norm of p is bounded uniformly in t. 
The L 2 -norms of 4>± can be estimated as 

/oo roo 
\cp±(t,u)\ 2 du = / \(p±(0,u)\ 2 du =: c, 
-00 J— 00 

where v! = u^ft. 

Near infinity, we work similarly with the fundamental solutions <f> ()3.2f) . Again using 
the Jost representation (see (pe2) and Lemma [3T3|) . we get terms depending only on tin 
as well as error terms which decay like 1/u and are thus in I?. 

The time derivatives can be treated in the same way, since a time derivative merely 
gives a factor of u which can be absorbed into h±. For the spatial derivatives we use 
similarly the estimates for the first derivatives of the Jost functions. ■ 



Our next step is to decompose the energy integral into a convenient form. To do this, 
we introduce a positive mollifier a E Cg°([— 1, 1]) with the properties a(— uj) = a(u) and 
J a{u))dw = 1. We define the function T(u — uj') by mollifying the Heaviside function 0, 

r(u-u/) = (e*a)(w-u/) . (8.3) 

We now substitute the Fourier representation of <J>h+ into the formula for the energy 
density (|8.2p . For simplicity we omit the indices no and H+ in what follows. Omitting 
the first positive summand in (|8.2|) . we get the inequality 

/OO /"OO 
dco / duj' e - i ( w - w ')* 
-00 J —00 
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sirr d A 



+T(u-u') (A d r ${Lu) d r ${u>) + sin 2 $ d cos v$(La) d cos #$(Lj')) (8.5) 
+(l_r)(w-a/) (A r $M dA(co') + sin 2 #d cos #$(u) d cos ^(cu^ j. (8.6) 



We multiply by a positive test function r)(u) S Cq°(R) and integrate over r and cos??. 
Integrating by parts in (|8.5p and (|8.6p to the right and left, respectively, we can use the 
wave equation 



d . d 1 ,, 2 ■-> , \i 



9 -,2 M 9 , 1 /„,.„ ; „2„o , ,x2 



■ sin 2 ?9— + (ao; sin 2 ■& + 



5 cos i? 3 cos 7? sin 2 $ 

to obtain 



= 



/CO /*1 / v f'OO f*l f'OO 

dr J dcos ■& rj(u) \J$3ft + dSBJ J = y dr y dcostf y efcj 



do/ «.-*(«"-"')* 



| -t/(u) (r 2 + a 2 ) (r a r $(o;') + (1 - T) d r $(u;) $(u/)) (8.7) 

+r?(n) (r^O + a-r)^)) 6(a;)i^y) 



where 

1 2 1 

g(u,r,<d) = — ((r 2 + a 2 )oj + ak) 5— (aw sin 2 # + fc) 2 , 

A sin d 

and we used that 

d ^2 _j_ q2 

— 7]( U ) = 7]'( U ) . 

We interchange the orders of integration of the spatial and frequency integrals and let 77 
tend to the constant function one. In the term corresponding to (|8.4j) . Lemma 18.31 allows 
us to pass to the limit. In ([8.71 I8.8j) the situation is a bit more involved due to the 
factors of T. However, since multiplication by T(u> — u>') corresponds to convolution with 
its Fourier transform T(t), we can again apply Lemma 18.31 and pass to the limit using 
Lebesgue's dominated convergence theorem. To make this method more precise, let us 
show in detail that the expression 

dr f dcostf f duo I drf e' 1 ^-^ 1 (\ - rj(u)j T(u - J) g{uj') ®(u) $(u') (8.9) 



'1 



tends to zero as rj converges to the constant function one. Rewriting the factor T with a 
time convolution, we obtain the expression 



/oo 
dr f (r) F(t) 
-co 



where F and f are defined by 



dr J dcostf (1 - 7](u)) $(t - r,r,i?) (g *W)(t - r,r,#) 



1 



CO 



f(r) = — / r(6) e~ ibT db . 
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Writing g as a polynomial in u, 



where 



go 



2u2 



a k 



k 2 



A sin 2 



2„Q ' 



gi = 2ak 



90 + gi u + 92 u 



3.10) 



r 2 + a 2 
A 



.92 



r 2 + a 2 ) 2 2 . 2 



A 



a 2 sin 2 tf, (8.11) 



the function g * <I> can be expressed explicitly in terms of $ and its time derivatives of 
order at most two. In order to compute T, we first note that the Fourier transform of the 
Heaviside function O is 

G(r) = L(-i?l+ v 6(T)), (8.12) 

where "PP" denotes the principal part. Using (18.3|) together with the fact that convolution 
in momentum space corresponds to multiplication in position space, we find that 



V(r) - ( -i ^ +tt5(t) j n(V) 



(8.13) 



where a is a Schwartz function with d(0) = (2ir) . According to Lemma 18.31 the func- 
tion F is uniformly bounded, 

\F(t)\ < sup\r]\ for all t e R. 

Using the rapid decay of T, we can for any given e choose a parameter L > such that 



/ 

J MS 



[-L,L] 



f (r) F(t) 



< e sup \r}\ 



On the interval [-L, L], on the other hand, the singularity of f at r = can be controlled 
by at most first derivatives of F, and thus for a suitable constant C = C(L), 



-L 



f(r)F(r) 



< 2C sup (|F| + \F'\ 

[-L,L] 



Using the rapid decay of $ and its time derivatives in u, locally uniformly in r, we can 
make supr^n (\F\ + \F'\) as small as we like. This shows that (|8.9|) really tends to zero 
as 7] goes to the constant function one. 
We conclude that 



£($£+) > / dr dcostf 



co roo 
did 

-oo J — oo 



dJ e -<(«-^)* 



k 2 a 2 k 2 



sin 



$(u;)$(u/) 



3.14) 



+ ^ 9 (u;') + (l-^) 5 H-min( 5 (a;), 5 (a; / )) < ^ ) H <^) ( w, ) • ( 8 - 15 ) 



We apply Lemma 18.21 to obtain 



(HHD > / dr / dcostf ( inf + 

7ri 7-1 V"^- 7 



fc 2 



a 2 fe 2 



sin 2 i? 



|*(i,r,0)|' 
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Using the explicit form of g, (|8.1(jp . we find that for sufficiently large J, 

mm > 4r / . dr / dcosti Mt,r,$)\ 2 . (8.16) 

2 J ri A J_ x 

It remains to control the term fl8. 15f) . We write the uj, cj'-integral of f|8.15j> in the form 

poo poo 

B := / duj / dw'e _i ( w - w ')*/i(w,w / )*H*(w / )> 



where 

h(u,u') = Fg(u') + (1-T)g(u) -mm(g(u),g(u}')) . 

Introducing the variables a = + oj') and 6 = |(u; — a/), and using that is a 

polynomial in w, a short calculation yields 

h(a + 6, a - b) = ( gi + 2g 2 a) 5(26) where 5(6) := 6 (e(6) - T(6)) . 

Using (|8.12| I8.13P together with the fact that the factor 6 corresponds to a derivative in 
position space, we obtain 

_ i an-2.a { r)\ = i f\ l(sT)ds 

K J 2vr dr \ r / dr J v y 

This is a smooth function which decays quadratically at infinity; in particular, it is inte- 
grable. 

We thus obtain for the Fourier transform of h the explicit formula 

-i roo poo 

h( T ,r') = 7A2 / dco I du'/i(w,w')e- i(uWt,) 



00 



= {gi5{T-r') + 2ig 2 8\T-r'))s(^^\ 
Using Plancherel for distributions, we obtain 

/oo /*oo 
dr / dr' ^(r,r') $(t - r) IfT^ 17 ) 
-00 J —00 

/oo _ 
dr5(r) 51 $(t - r) $(t - r) 
-00 



+*SSJ - r) $(t + t) - $(t - t) d t $(t + r) 

Integrating over space, we can use the explicit formulas for g\ and g 2 and apply Lemma [8.3l 
to obtain 

poo pi poo 

(gI5D = dr dcostiB > -C(n ,*o) / |5(r)| dr . 

Jri ■/— 1 J —00 

We now let a tend to the Dirac delta; then a tends to the constant function (2-7r) _1 . As 
a consequence, the L 1 -norm of 5 tends to zero, and thus (|8.15p becomes positive in this 
limit. Hence the energy is bounded from below by (I8,16h , This concludes the proof of 
Theorem 18.11 
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8.2 Pointwise Estimates for the Low-Energy Contribution 

The low-energy contribution can be written as 

2 



-V 

2vr ^ 



duj 



00 



n=r -uu 0,6=1 

We now derive pointwise estimates for the large angular momentum modes. 

Theorem 8.4 For any uq < u\ there is a constant C > such that for all oj £ (—2 J, 2 J) \ 
{loq,0} and for allu,u' G (uc^i), 



£ 

71=1 



2 ^fM <^V) 



a,6=l 



< c. 



.17) 



Proof. From (|5.5p the coefficients t a & have the explicit form 



1 + Re 



(t 



a 
P 



Tm 



ab 



where the transmission coefficients a and /3 are defined by 

<fi = a cj> + (3 <fi , 

and (pi = Hecf), 4>2 = Im0. The estimates of Section [4.31 are obviously valid for u> in any 
bounded set; in particular for u> G (— 2J, 2J) \ {wo,0}. We use these estimates in what 
follows, also using the same notation. We choose n\ so large that u + < uq, and thus on 
the whole interval (uq, 2ui) the invariant disk estimates of Lemmas 14.21 and 14.81 hold. 

Rewriting the expression t a \)(j) a (j) h with the Green's function (see the proof of Lemma fS. 1 1) . 
this expression is clearly invariant under the phase transformation cf) — ► e <f>. Thus we can 
arrange that 4>(2ui) is real. Then the transmission coefficients are computed at u = 2u\ 

by 



We thus obtain 



a 



2 liny 



ji=2iji 



/9 



Hence 



1 + 



limy I 

2 7 ^7 < 



2 <^>Imy 
4|0| 



(y - y) 



u=2u\ 



\y - V\ P 2 Re (y - y) ' 
where all functions are evaluated at u = 2u\, and where we used the relation 

p = J — - ' 

Imy 



3.18) 
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which is an immediate consequence of (|4.9p and w((f>, 4>) = 2iQ. From Lemma l4T71 and (|4.42j) 
we know that for A sufficiently large, 

1 



Re (y - y) > 



C 



Using the above formulas for t a b, we conclude that 



(a,6)^(2,2) 



6 v: 



< 12 C 



p(u) p{u') 



p{2u x ) /o(2ui)' 



The argument after (|4.43|) shows that the two factors on the right decay like exp(— vA/e). 
It remains to consider the case a = b = 2. Taking the imaginary part of the identity 



4>{u) = (f>(2ui) exp 

and using that 4>{2ui) is real, we find that 

p(u) = p{2u\) exp 

and thus 

1 02 0)| = p{u) 



2 a i 



Re y 



2ui 







/•2ni 


sin Imyj 


< P{u) 


/ Imy 






J u 



From (|4.40p we see that Rey is positive on [u, 2u\]. Using the relation p'/p = Rey, we 
conclude that p is monotone increasing, and (|8.18p yields that Imy is decreasing. Hence, 
again using (|8.18|) . 



|02O)| < p(u) Imy(u) (2ui - u ) < Imy(u ) {2u\ - u ) 

Using the above estimates for a/ (3, we conclude that 

1 



n 



t22 P{u)<j> l (v!) 



< 3Imy(u ) (2u\ - u ) . 



(8.19) 



The invariant region estimate of Lemmas 14.21 and 14.81 yield that 

7 fUo 



Imy(ito) < c\Q\ exp 



and we conclude that (|8.19p again decays like exp(— y/X/c). 

Since the eigenvalues A n scale quadratically in n (|2.1ip . the summands in (|8.17p decay 
exponentially in n. ■ 



This theorem gives us pointwise control of the low-energy contribution, uniformly in 
time and in no, locally uniformly in space. To see this, we estimate the integral represen- 
tation KHZ} by 

— tT$ a (r ■&) <^ b *n> 



r,0)| < -J2 / *»\Xl{u)\ E 

Z7r n=1 J-oo b=1 



In order to control the factor to 1 , we write the energy scalar product on the right in the 
form (2.15)], which involves an overall factor u. Now we can apply Theorem 18.41 
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8.3 Decay in L 



In this section we complete the proof of Theorem I l.li Let K C (ri, oo) x S 2 be a compact 
set. The L 2 -norm of *™° can be estimated by 

\\<S> n °(t)\\LHK) < \\*Z(t)\\l*(K) + \\K°-(t)\\LHK) + m°(t)\\L H K) ■ 

According to Theorems 18.11 and 18.41 these norms are bounded uniformly in no and t. 
Furthermore, our estimates imply that the sequence * n °(i) forms a Cauchy sequence 
in L 2 (K). To see this, we note that for any ni,ri2, 

- * n2 || L2( x) < - *S 2 + ) + - *£ 2 _) + II*? 1 - ^ 2 \\LHK) , 

uniformly in t. The energy terms on the right are the sums of the energies of the individual 
angular momentum modes. All the summands are positive due to Theorem 18. 1\ and thus 
the energy terms become small as n\ y n2 — * oo. The same is true for the last summand 
due to our OD estimates of Theorem 18.41 We conclude that * n °(i) converges in L 2 oc 
as no — > oo, and the limit coincides with the weak limit, which from [5] we know to be the 
solution *(t) of the Cauchy problem. 

To prove decay, given any e > we choose no such that ||*(t) — * no (i)||L2(A-) < £ for 
all t. Since * n ° is continuous in u and has rapid decay, uniformly on K, the Riemann- 
Lebesgue lemma yields that * n °(£) decays in L°°(K) C L 2 (K). Since e is arbitrary, we 
conclude that *(i) decays in L 2 (K). 

Applying the same argument to the initial data H n ^Q, we conclude that the partial 
derivatives of *(t) also decay in L 2 {K). The Sobolev embedding H 2 ' 2 (K) ^ L°°(K) 
completes the proof. 
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